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Abstract. The hidden super symmetry of the monopole found by De Jonghe et al. is 
generalized to a spin \ particle in the combined field of a Dirac monopole plus a A 2 /r 2 
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potential [considered before by D'Hoker and Vinet], and related to the operator introduced 
by Biedenharn a long time ago in solving the Dirac- Coulomb problem. Explicit solutions 
are obtained by diagonalizing the Biedenharn operator. 
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1. Introduction 

In [1], De Jonghe et al. found that a spin-| particle of mass, here chosen m = h, 
and charge q in the field of a Dirac monopole of unit strength, described by the Pauli 
Hamiltonian 

|: d-i) h d = 

admitted a 'hidden' supercharge Qd given by 
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(1.2) Qd = (a ■ £+1), 

where £ is the [non- conserved] orbital part of the angular momentum, £ = rx7?, n = p—qA. 



Their discussion is based on the study of Killing- Yano tensors. The supercharge Qd 
anticommutes with the supercharge Qd found by D'Hoker and Vinet [2], where 

(1.3) Q D = o ■ 7T. 

These supercharges form, together with the [conserved] total angular momentum 
J = L + ^, L = £ — qr/r, a closed, non-linear algebra [1]. 

In this Letter, we (i) generalize the result of De Jonghe et al.; (ii) relate it to earlier 
work of Dirac [3], Biedenharn [4], and Berrondo and Mcintosh [5]; (iii) use it to solve the 
system, and (iv) discuss a full Minkowski space generalisation. 

Very recently [6] , Plyushchay discussed related problems, but from a rather different 
veiwpoint : while our results here are derived from supersymmetric quantum mechanics, 
he uses pseudoclassical mechanics with anticommuting (Grassmann) variables. See also 
the comments on [1] in [7]. 
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2. The generalized monopole system 

Let Ad denote the vector potential of the Dirac monopole of unit strength. Let A be 
a positive constant, q > a half-integer, and consider, setting = A/r and A = qAo, 
the gauge field A a , (a = 1, . . . , 4) on R 4 . We use the Dirac matrices 

< 2i » ?=U "*)• '=U l2 )- 5 =( 12 -i.)- 

When restricted to fields which do not depend on x 4 , the associated Dirac operator 
ft = 7 M .D M , Dp = — iAp, reads 

")-(....„ '-*)■ 

^ anticommutes with the chirality operator 7 5 . 
The square of the Dirac operator, 

is a supersymmetric Hamiltonian [8]: the non-zero-energy parts of the chiral sectors, 
defined as the eigensectors of 7 5 , are intertwined by the unitary transformations 

(2.4) U = T— = and U' 1 = U ] = — = T t , UH X U ] = H . 

yHi y/Hi 

The partner hamiltonians Hi and Hq, have therefore the same positive spectra, and the 
spectrum of ft can be obtained from that of ft 2 . For further analysis, it is convenient 
to set 

(2.5) o r = a- (f/r), z = a-£+l. 

Note that a 2 = 1, and that z anticommutes with o> and a • tt, {z, a r } = and 
{z, a ■ n} = 0. Also z is equal to the supercharge Qd of (1.2). For A ^ 0, z satisfies the 
relations zT = —T^z, zT^ = —Tz, zH\ = H z and zHq = H\z. Therefore, the operator 

commutes with the Dirac operator ft , and hence also with its square. Using (a ■ L) 2 = 
L 2 + id ■ (L x L) = L 2 — a ■ L, one proves furthermore that K? = z 2 = J 2 + 1/4 — q 2 . 
Thus, since the eigenvalues of J 2 are j(j + 1), j = q — 1/2, q + 1/2, . . ., the operators z 
(and /C) have irrational eigenvalues, 



(2.7) 



« = ±V(j + V2) 2 - 4 2 



Monopolesusy 



3 



The operator /C is hermitian because j > q — 1/2. For the lowest allowed value 
j = q — 1/2, however, k vanishes and /C is not invertible. The operator /C has been used 
by Dirac in the study of the relativistic hydrogen atom [3] long time ago; its form adapted 
to the monopole, (2.7), was found by Berrondo and Mcintosh [5]. It is more convenient 
to use, however, the hermitian Biedenharn operator [4] 

(2 ' 8) T --{ a-£ + l-\o- r 

Here y = z + \a r and x = z — Acr r , so that xa r = —a r y. The eigenvalues of T, 




(2.9) 7 = ±y/(j + 1/2) 2 + A 2 - q 2 , (sign 7 = -sign k) 

are in general still irrational. However, owing to the presence of A 2 , the operator F is 
invertible whenever A 2 > 0. The clue is that, in terms of V, if) 2 becomes simply 



(2.10) 



1x2 , r(r + i) 



Here p r = — i(l/r)d r r = —i(d r + 1/r) is the hermitian operator conjugate to r. This is 
conveniently checked by writing, using the radial form a • tt = —ia r (d r + 1/r — z/r). The 
supercharges T and as 

T = -lay. (d r + i - ^) = -i (a r + i + f )a r , 
(2.11) , 

Tt = -za r ^ r + i-f) = -z(a r + 1 + ^)a r . 

The self-adjointness of p 2 requires |A| > 3/2 [8]. It follows from our previous 
formulas that V anticommutes with if) and commutes therefore with its square. Hence the 
shifted operators x and y commute with the partner hamiltonians, [x, H ] = = [y, Hi]. 
In conclusion, our 4-component operators satisfy the non-linear algebra 

P 2 = H, {r,^}=0, [T,H]=0, T 2 = f + l/4 + \ 2 -q 2 , 

[J, r] =0, [J,p]= 0, [J, H] = 0, [Ji, J,] = ie ljk J k , 

which only differs from that found for the 2-components objects of [1] in the appearance 
of A 2 in T 2 . Note that it is now if) which plays the role of Qd, and V plays that of Qd- 

3. Explicit solutions 

A look at (2.10) shows that the Biedenharn operator V plays clearly a role analogous 
to angular momentum. Since T and J commute, they can be simultaneously diagonalized. 
A convenient basis is found as follows [9]. Let us first assume that j > q + 1/2, and let 
L± = j ± 1/2 be the orbital angular momentum quantum number. Consider first the 
two-component spinors 



<3,) w^^i: 
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where the Vs are the monopole harmonics, defined in [10]. The <^'s satisfy 

J 2 = j(j + 1), J3 = fj,, (/i = —j, • • ■ , j), L 2 = L±(L± + 1), and the action of o r upon the 

</?'s can be obtained from [11]. Then, dropping the upper index n, and using 

a ■ £ = J 2 — L 2 — 3 / '4 + qa ■ r/r , it may be shown that the 2-spinors 



(3.2) 



X± = 



q 



j + 1/2+ I K 



satisfy z x± = ±|/c| X± and a r x± = Xt> as wel1 as ^ = JU + 1 ), J s = V, = ~h ' ' -J)- 
Hence, defining 7 by 7 2 = k 2 + A 2 , we have 



(3.3) 



<t>A 



l«l + M) x+ - A x-> 0- = A x+ + (|«| + ItDx- 



= (l«l + ItI) x+ + A x_, = -A x + + (|«| + ItDx- 



diagonalize a; and y, x0± = =F I 7 I 4>± an d 2/$± = -F I 7 I The operator a r = a ■ r/r 
interchanges the x and y eigenspinors, a r <j>± = $ T , a result which also follows directly 
from xa r = —o r y. 

Whenj = q—\, there are no L_ = q— 1 states, though: we only have 2(g— 1) + 1 = 2g 
states with L = L + = q , namely 



(3.4) 



Thus, for j = q — \i no is available, and x+ = —X- = is annihilated by z. 
Therefore, there are no 0_-states in the 7 5 = —1 sector, and no $ + states in the 7 5 = 1 
sector. In each 7 5 sector, (3.4) yields in turn (2q) states, namely 



(3.5) 



They are eigenvectors of x and y with eigenvalues ±A, respectively, and are still inter- 
changed by a r . In conclusion, the eigenfunctions of if) 2 are 



u± 



u± 



(3.6) 



// 



where u± solves the radial equation 









</>+ 



for 




for j > q + 1/2, 



for 



T 



7 5 = -1 



for j = q — 1/2. 



(3.7) 



IV , 7(7 + 1) 



a r + - + 



E 



u±(r) = 0. 
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There are no bound states; the scattering states involve Bessel functions: 
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(3.8) u ± (r) = -Lj h+h (VEr). 



4. Symmetries 

A spin particle in the field of a Dirac monopole has an o(2, 1) symmetry, generated 
by H = Hq = 7? 2 and by dilations and expansions [12], 

(4.1) D = tH- 1/2{tt -r) K = -\t 2 H + tD + \r 2 . 

This symmetry has been extended to the Pauli Hamiltonian (1.1) with formally the same 
generators (4.1), with H = Hd replacing H [2]. The supercharge Qd is a square-root of 
Hd- Commuting Qd with the expansion, K, yields a new fermionic generator, namely 

(4.2) S = i[Q D , K] = -La-(r- ift), 

v2 

and it is then readily proved that the bosonic operators H D , D, K close, with Q and S 
into an osp(l/l) super algebra. [2, 8, 13]. Now remarkably 

(4.3) i[Q D ,S\-± = z, 

and z 2 is a Casimir operator of this osp(l/l) [8]. 

The same bosonic o(2, 1) symmetry arises for the generalized monopole system (2.4). 
The Dirac operator Q = if) is a square-root of H by construction. However, 



(4.4) Q* = 7 5 Q = 



-> - -A 



is a new square-root, {Q*, Q*} = H. Commuting K with Q and with Q* yields 

S = 7 5 7 ■ f - tQ and S* = -rfS. 

In this way, we get two, independent, super-extensions of the bosonic o(2,l). The two 
osp(l/l)'s do not close yet: the "mixed" anticommutators between the Q-type and S-type 
charges yield a new bosonic charge, namely 

(4.5) Y = {Q, S*} = -{Q\ S} = 7 5 (z + \) - Xa r , 



that commutes with the other bosonic charges. The four operators H,D,K,Y do close 
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finally with the four fermionic charges Q,Q*, S, S*, 



(4.6) 



[Q,D] 


= iQ, 


[Q*,D] 


= iQ\ 


[Q,K] 


= -iS, 


[Q\ K] 


= -iS* 


[Q,H] 


= o, 


[Q\H] 


= o, 


[Q,Y] 


~iQ\ 


[Q\Y] 


= iQ; 


[S,D] 


= -iS, 


[S*,D] 


= -iS* 


[S,K] 


= o, 


[S*,K] 


= o, 


[S,H] 


= 2zQ, 


[S*, H) 


= 2iQ\ 


[S,Y] 


= -iS*, 


[S*,Y] 


= iS, 


{Q,Q} 


= H, 


{Q*,Q*} 


= H, 


{S, S} 


= 2K, 


{S*,S*} 


= 2K, 


{Q,Q*} 


= o, 


{S, S*} 


= o, 


{Q,S} 


= -D, 


{Q*,s*} 


= -D, 


{Q,s*} 


= Y, 


{Q\ s} 


= -Y. 



which are the commutation relations of the osp(l/2) superalgebra, to which spin adds an 
extra o(3) [8]. Now the Casimir of osp(l/2) is the square of 



(4.7) 



i[Q,S\-i=i[Q*,S*]-i 



which provides a nice interpretation for the Biedenharn operator V. Similar algebras were 
studied in [14]. 

5. Particular cases 

(i) For A = , we have Q = = Qd-, Hi = Hq = Hp, the Pauli Hamiltonian in 
a pure monopole field [2]. The 4-component Hamiltonian is simply diag(i?£), Hp); the 
Biedenharn and the Dirac operators are related as V = —i-f 4 JC. In this case, we recover 
the formulas in [1, 2] . 

(ii) Another particular value is \ = ±q, when the situation is similar to that in Taub- 
NUT space [15]: the spin drops out in one of the chiral sectors, while the Pauli term gets 
doubled in the other. For A = q, for example, the Hamiltonian (2.4) reduces to 



(5.1) 



H 



Hn 



-2q- 



a ■ r 



zr 



- 2 + q -i 



Here Hq describes a spin particle in the combined field of a Dirac monopole and of an 
inverse- square potential, while H\ corresponds to a particle with anomalous gyromagnetic 
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ratio 4. The Biedenharn operator has half-integer eigenvalues, 7 = ±(j + \). Note that 
the 7 5 = —1 eigenspinors now reduce to those in Eq. (3.1), <fi± oc ip±. 

Assume first that j > q + 1/2. Since 7(7 + 1) is now the same for — I7I as for I7I — 1, 
these values lead to identical solutions. Therefore, in each 7 5 sector, the energy levels are 
two-fold degenerate. The numerator of the r~ 2 term reads 

( , 2) « + M + *> fOT ( 7> ° 

so that we indeed get the same equation with 7 > for j, as with 7 < for (j — 1), 
provided that (j — 1) states do exist. In both cases, (3.8) yields r" 1 / 2 J j+1 (y/Er). 

The two-fold degeneracy is hence also explained by an extra o(3) symmetry in addi- 
tion to the rotational symmetry: spin is trivially conserved for H , and this is exported to 
Hi by supersymmetry. The extra o(3) symmetry is generated hence by the spin vectors 

(5.3) £ = \o for H , Si = S U for H u 

where U and U -1 = are the intertwiners of (2.4). The two-fold degeneracy corresponds 
precisely to this o(3) symmetry. For j = q — 1/2 , half of the states are missing. 

The system admits further symmetries. Firstly, Hq admits the non-relativistic con- 
formal o(2,l) symmetry [10]; supersymmetry exports this to the partner Hamiltonian 
Hi. The symmetries combine with if) and —i^lf) into an osp(2, 1) superalgebra [8] . 

(iii) Replacing the scalar potential A/r by q(l — 1/r) — which corresponds to the 
long-range limit of the scalar field of a self-dual "BPS" monopole [14] — we get 



(5.4) 



#1 

Hn 



/ ^ 9 2q 2 q 2 9 n a ■ f \ 
1 7? 2 - -2- + \ + q 2 - 2q 



9 ' 1 "1 o 



2 2q 2 q 2 2 



\ f T / 



The "lower" Hamiltonian, Hq, has again gyromagnetic ratio yielding an extra o(3) 
symmetry. Its properties have been thouroughly studied by Mcintosh and Cisneros, and 
by Zwanziger [17], who found that it admits a Kepler- type dynamical symmetry. Its 
superpartner Hi describes a spinning particle of anomalous gyromagnetic ratio 4: this is 
the 'dyon' of D'Hoker and Vinet [18]. Then supersymmetry can be used to transfer the 
symmetries of H to Hi [16]; the system can be solved using the Biedenharn method [9]. 

6. Minkowski space extension 

In this section we turn to the Dirac-equation in Minkowski space-time, with a com- 
bined Coulomb-monopole and massless scalar background. This relativistic system of 
equations has properties very similar to the Euclidean Dirac problem considered above, 
and can be solved by analogous methods. 

In Minkowski space-time, with metric i]^ = diag(— 1, +1, +1, +1), we use Dirac 
matrices with the properties {7^ , 7^} = 2^ and 7 M t = -707^70- 
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Our starting point is the Dirac equation 

(ift + m + gip)ip = . 



(6.1) 



Here ip = —g/r is a dynamical scalar background and m the mass (which can be taken 
as the vacuum expectation value of the scalar field). Eq. (6.1) describes the motion in 
the long-range field of a Julia- Zee dyon [19]. The covariant derivatives contain electro- 
magnetic background potentials corresponding to a Coulomb field 

D = d + iq(p , 4> = q/r , 

and the field of a magnetic monopole taken, as in Sec. 2, to be of unit strength. 

Multiplication by (—ift + m + gip) gives a Klein-Gordon type equation, which on 
stationary states ip(r,t) = ex.p(—iEt)ipE(r) takes the form 



— (E — qcj)) 2 — (V — iqA) 2 + (m + gp) 2 - iqa^F^ - i 91 ^ <p 



E 



0. 



(6.2) 



The generalized Klein-Gordon operator can be block-diagonalized in (2 x 2) blocks, 
with 2-component eigenspinors ip± satisfying 



-(E - qcj)) 2 - (V - iqA) 2 + (m + gip) 2 - a ■ (qB ± VA) ^ E ,± = 0, (6.3) 



with A = \/g 2 (p 2 — q 2 (j) 2 . Note that the square root is real for g 2 (p 2 > q 2 (p 2 , and 
imaginary for g 2 (p 2 < q 2 (p 2 . For the case of the Coulomb and scalar potentials this 
becomes 

A = - , A = vW 7 ??. (6-4) 



Defining the operators I and J as in Sec. 2, we may cast (6.3) into the form 



1 



l2 



£. a-q 2 T Aa r ) — (E — ^) 2 + (m - —) 2 ]^e,± = 0. (6.5 



To make contact with the Biedenharn operator and the work of Sec. 2, we note that the 
operators y and x of (2.8) occur as blocks in (6.5). Writing here r + = — y and T_ = — x, 
we have 

r± (r± + 1) = P + A 2 - (q ± A) o T , (6.6) 



and can hence rewrite (6.5) in the standard form 

2 



- d r + 



1 



+ ^r ± (r ± + i) + 



2k 



+ e 



E,± 



m 



(6.7) 



E 2 . Note the 



Here the constants k and e are given by k = qqE — ggm and e 
symmetry under the simultaneous exchange of (E, qq) <-> i(m, gg). The eigenvalues of r± 
are the ones given in (2.10): here 



7 = ±^(j + 1 1 ) 2 + g 2 ~g 2 -q 2 (q 2 + l) 



(6.8) 
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Introducing the notation Z 7 , where Z 7 = 7 for 7 > and Z 7 = — (1 + 7) for 7 < 0, we see 
that the eigenstates of T± satisfy the equation 



1\ 2 1 2k 



rl> Em = 0. (6.9) 



The spectrum of eigenvalues for bound states is well-known from atomic physics: 

e = k 2 /n 2 f , n 7 = l-M 7 + iV, N = 0,1, 2, ... . (6.10) 
However, in this case the bound-state energy eigenvalues themselves then are given by 

E(j,N) = m (ggqq ±n lS Jn 2 + q 2 q 2 - g 2 g 2 ^j . (6.11) 

For the ground state j = 0, N = 0, the wave equation factorizes as expected on the 
basis of supersymmetric quantum mechanics: 



-d r + ^-^ + -^- a r + -^-^ + ^ Vo = , (6.12) 
r la) \ r la) 

with 7 fl and k g the ground state values of 7 and k. 
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